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Local helioseismology and the active Sun 
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The goal of local helioseismology is to elicit three-dimensional information about the sub-surface (or far-side) structure 
and dynamics of the Sun from observations of the helioseismic wave field at the surface. The physical quantities of interest 
include flows, sound-speed deviations and magnetic fields. However, strong surface magnetic fields induce large perturba- 
tions to the waves making inversions difficult to interpret. The purpose of this paper is to outline the methods of analysis 
used in local helioseismology, review discoveries associated with the magnetic Sun made using local helioseismology 
from the past three years, and highlight the efforts towards imaging the interior in the presence of strong magnetic fields. 
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1 Introduction 

Helio- and astero-seismology provide the means to "see" 
beneath the surface of the Sun or star, and to gain an un- 
derstanding of the physics of stellar interiors. The analysis 
of the frequencies of the global modes of oscillation allows 
us to study the internal stellar structure only as a function 
of radius and unsigned latitude. With high spatial resolution 
observations of the velocity signal on the Sun's surface, lo- 
cal helioseismology is able to image the three dimensional 
subsurface structure of the Sun and image magnetic activity 
on the far-side. Local helioseismology aims to interpret the 
full wave field observed at the solar surface, not only the 
frequencies of the normal modes. We restrict our attention, 
here, to local helioseismology. 

The Doppler signal, which is the line-of-sight velocity 
of the solar surface, is the fundamental observable. Full-disk 
high-quality data is obtained from instruments including the 
Michelson Doppler Image (MDI) onboard the Solar and He- 
liospheric Observatory (SOHO) (Scherreret al. 1995), the 



Globa l Oscillation Network Group (GONG) dHarvev et al. 



Il996i). and the Helioseismic and Magnetic Imager (HMI) 
onboard the recently launched Solar Dynamics Observa- 
tory (SDO). MDI also runs higher-resolution campaigns on 
smaller fields-of-view. The Magneto-Optical filter at Two 
Heights (Fijisterle_et al. 200 4) observes simultaneously at 
two heights in the atmosphere allowing an analysis of the 
vertical propagation of waves. 

Various local helioseismic techniques exist to analyse 
the oscillations. Th e main ones are time-di stance ana lysis 



mode helioseismology. By analysing frequency shifts over 
small regions of the solar surface, the direction and ampli- 
tude of subsurface flows can be determined. Time-distance 
analysis computes the travel time of a wave packet trav- 
elling between two points on the solar surface. Helioseis- 
mic holography (acoustic imaging) reconstructs the subsur- 
face wave fields by propagating waves either forward or 
backward in time. Other methods include Fourier-Hankel 



analysis (IBraun et al.lll987h . desi gned sp ecifically to anal- 
yse waves surrounding a sunspot. IWoodard (.2009.) used di- 
rect modelling to interpret the wave correlations in wavevector-| 
frequency space. 

The forward problem is then to compute the expected 
seismic observations from a particular model of the Sun. 
The inverse problem is to infer the internal properties of the 
Sun from the seismic observables. Local helioseismology is 
the only way to elicit the three-dimensional structure of the 
solar interior, including structures and flows on a range of 
scales. 

In this paper we focus on the local helioseismology of 
sunspots, active regions and large scale flows varying with 
the solar cycle. There are particular problems associated 
with doing helioseismology in the immediate vicinity of 
strong surface magnetic fields. The forward problem used 
to infer the subsurface properties has previously been based 
exclusively on solving hydrodynamic equations in the ab- 
sence of a magnetic field. Additionally, the actual inversion 
method assumes that any perturbation to the waves is small, 
however it has been demonstrated that strong surface mag- 
netic fields cause large perturbations to the waves. Know- 
ing the subsurface structure of sunspots and active regions 



(IDuvall et al.lll993h . ring diagr am analysi s (Hill 1988) and 
acoustic imaging/holography (.Chang et al..l997:.Lindse v & BrMf ""^^"^i" s^"sP°t "^^'^els ^"^^ hopefully determine their 
19971). Ring-diagram analysis is most closely related to global deepei" ^^""'^^^i"^' ^he immediate need now is to quantita- 

tively model wave propagation through magnetic fields and, 
subsequently, account for the effects when inverting for the 
subsurface structure. 
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An important question in solar physics is to understand 
the mechanics of the solar dynamo. Local helioseismology 
has identified large scale flow variations that appear to be 
intimately connected to the surface activity. Understanding 
such connections will help to decipher the larger scale dy- 
namics of the magnetic Sun. Other goals of local helioseis- 
mology are to measure the meridional flow and changes at 
the base of the convection zone, which would help constrain 
models of the dynamo. 

This short review will cover the basic techniques used in 
local helioseismology in Section |2] How the interior struc- 
ture is subsequently inferred is described in Section[3] Sec- 
tion|4]will then cover some of the results, over the past three 
years, from local helioseismic analysis tied to solar activity. 
A concentrated effort towards modelling waves propagating 
through magnetic fields to understand the observations will 
be reviewed in SectionB) Rec e nt reviews of local hehoseis - 
mology inclu de BirchI (l2008h : iGizon & ThompsonI (l2007b : 
Komm et al.l (1200 6). For more e xtensive reviews of local he- 
lioseismology see, for example. iGizon et al.l (120101) . 



and a longer travel time if the flow is in the opposite di- 
rection. One way to extract the travel times from the cross- 
covariance is to fit Gabor wavelets to the cross-covariance 



dPuvall et al.l Il997h as a function of time-lag. The cross- 
covariance ca n also be fit using a one parameter fit lGizon & Birch[ 
( I2OO2I I2OO4I) . Often, the cross-covariance is computed be- 
tween a central point and a surrounding annulus, and aver- 
aged over the annulus. This technique is used to study the 
divergence of flows. The direction of th e flows is gained b y 
dividing the annulus into quadrants (e.g. Duvall et al.ll996h . 
Filtering the observations is a common practice. In ridge 
filtering only modes of a particular radial order are retained, 
and in phase-speed filtering waves within a certain range 
of wave speed (or depth penetration) are retained. It has re- 
cently come to light that the filtering process has a strong ef- 
fect on the travel-time measurements and must be included 
in any interpretation of the observations (.Braun & BirchL2008l :[ 
Gizon et alJl2009HMoradi et al.ll2009bl) . 



2 Analysis techniques 

2.1 Ring-diagram analysis 

Typically ring-diagram analysis uses many small regions 
covering a significant fraction of the visible solar surface. 
Beneath each patch, two independent quantities, e.g. sound 
speed and density, may be inferred as a function of depth 
from which other quantities can be calculated. By combin- 
ing the information in each patch, a three-dimensional pic- 
ture is achieved. Typically this method is used to probe the 
top 30 Mm of the convection zone. 

For each patch a three-dimensional power spectrum is 
computed. A cut of the power spectrum at constant fre- 
quency results in an image of rings (.HilL1988:) . 

Variations in the power spectra indicate particular sub- 
surface properties. The presence of a flow shifts introduces 
a Doppler frequency shift and affects the shape of the rings. 
The degree to which the different modes are affected gives 
some indication of the depth and strength of the flow. Cuts 
at constant wavenumber can also be considered, where the 
power spectrum is a function of azimuthal angle and fre- 
quency. The measurement technique used in ring- diagram 
analysis and their interpr etation are discussed in (e.g. lBasu et al 
I999l:lHaberetaDl2000l) . 



2.3 Acoustic imaging/holography 

Acoustic imaging ("Chang et al.l Il997h and helioseismic 
holography ( Lindsev & Braun 119971) take the surface 
Doppler observations and reconstruct the subsurface wave 
field from ingoing and outgoing waves with respect to a 
pupil. The ingoing waves are propagated forward in time 
and the outgoing waves are propagated backward in time 
using Green's functions and a knowledge of the solar inte- 
rior Comparisons between the forward and backward prop- 
agating waves are used to derive pha se shifts and travel time 
variations. Far-side helioseismology (ILindsey & Braunl2000l) [ 
was developed using this method, although it is also possi- 
ble with time-distance techniques (IZhadl2007h . 



2.4 Fourier-Hankel 

Specifically designed with the intention to s tudy waves sur- 
rounding sunspots, Fourier-Hankel analysis (IBraun et alJl987h [ 
decomposes the observed oscillations into in- and out-going 
waves in a cylindrical coordinate system centred on a sunspot.| 
The observed waves are observed in an annulus surrounding 
the sunspot. From this the in- and out-going waves can be 
compared, in particular, their amplitude and phases dBraun 
I995l:lBraunetaDll988h . 



Inversions 



2.2 Time-distance 



Time-distance helioseismology (IDuvall et al.ll993h involves^ 
measuring the (phase) travel-time between two points on 
the solar surface from the cross-covariance function of the 
Doppler signal at the two points. The cross-covariance is 
sensitive to the wave speed in the solar interior. For exam- 
ple, if there is a horizontal flow in the direction from point 
1 to point 2 then this will result in a reduced travel time. 



The inverse problem in helioseismology involves extracting 
quantitative physical information regarding Solar subsur- 
face inhomogeneities; such as sound-speed, density, mass 
flows, and magnetic field. This requires accurate and realis- 
tic forward modelling of the solar osciUations, given a pre- 
scribed background mod el. Model S 

dChristensen-Dalsgaard et al.l Il996h is often used for this 
purpose in helioseismology. Furthermore, robust and reli- 
able helioseismic travel-time measurements must be com- 
bined, or averaged, in some sensible way in order to best 
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exploit the physical quantities of interest. It is also very im- 
portant that the noise propertie s of the observations be well 
understood and accounted for (IGizon & Birchll2004l) . 

Heretofore, the effect of the inhomogeneities on the os- 
cillations has been assumed to be weak. This ultimately 
yields a linear integral relationship relating observations, 
3D sensitivity kernel functions, and the subsurface perturba- 
tions. As yet, no non-linear inversions have been attempted. 

There exist two primary procedures for computing in- 
versions for time-distance and ring-diagram a nalysis; The 
Regu larised Least Squares (RLS) method (e.g. iKosovichev 



1996r) and the Optimally Local i sed Averaging (OLA) 



It is known that active regions are surrounded at the sur- 
face by a general inflow and when averaged longitudinally 
these flo ws introduce a perturbation to the meridional cir- 
culation. I Spruit! (120031) suggested that the inflow is due to 
the increased radiative loss of the active region compared 
to the surrounding Sun, and this can explain the solar cycle 
variations of the meridional flow. In agreement with this, it 
has been shown bv lOizon & Rempell (2008.) that a model in- 
cluding increased radiative loss in the active regions can re- 
produce the inflows and explain the observed surface solar- 
cycle variations. 

B y mask i ng ou t active regions and their surrounding 



method (e.g. jj^kiewicz et al.l20()7tlPiipers & Thompsonll992l^ w^'lgiZOI]| (l2004|) finds that the influence on longitudinal 

. • • J J 1 -1 1 r averagesof active region flows is significantly reduced. This 

Accurate mversions depend heavily on the accuracy of ° ° o j 

the kernel functions, which must reflect all the aspects of 
the measurement procedure, including filtering. The ray ap- 
proximation can be used to calculate the kernels 
(IKosovichev 1996) . but it is only useful if the background 
properties do not vary on scales smaller than, or compara- 
ble to, the wavelength. The first-order Born approximation 



implies that the flows surrounding active regions are the ma- 
jor contributors to th e time varying component of the merid- 
ional flow. Recently. lGonzalez Hernandez et aL (.2010.) have 
shown that in this extended solar minimum the residual merid-l 
ional flow developed before the onset of visible magnetic 
activity, suggesting that there is still a component indepen- 
dent of the presence of active regions. 



offers a more reliable way to compu te kernels, since it ac- 
co unts for finite wavelength effects ( Birch & Felder 20041 

fzon & Bi rchl2"002l:lHung et al.l200lb . Recentlv.l Birch eTall 4.3 Torsional oscillations 
m have computed sens itivity kernels for ring-diagram ^^^^.^^^j ^^^^^^^.^^^^ ^^ ^ands of sHghtly faster rotation 
analysis and|Birch&Gizon|(l2007l) have computed Bornker-I ^^i^^ijing with the activity belt which migrate towards the 

nels for time-distance hehoseismology. . , i t-, r, , , ■ t a ^u^u 

^■' equator over a solar cycle. The flow band associated with the 

new solar cycle has been found to be moving more slowly 

than in the lead up to the previous cycle (Howe et al. 20()9|), 

possibly providing a clue as to why the onset of Solar Cycle 

24 is later than expected. 



Jackiewicz et al.l ( 120081) have recently obtained maps of 



supergranular vector flows in the near-surface layers using 
an OLA inversion of travel times 



4 Recent observational results 

4.1 Meridional flow 

Meridional flow is the flow from the Equator towards the 
poles near the surface, with a return flow at some depth. The 
meridional flow is important for magnetic flux transport in 
solar dynamo models. An estimate of the meridional flow at 
the base of the convection zone predicts an amplitude of a 
few m s~^ whereas the sound speed at this depth is a few 
thousand m s~^, maki ng it difficult to detect as discussed by 



Braun & BirchI (|2008|) 
r 



■n 



Recently Gizon & Rempel (120081) and 

[Gonzalez Hernandez et al. ( 2008) have measured the mean 
amplitude of the flow using time-distance and ring-diagram 
techniques respectively and get consistent values of 10 — 15 
m s— 1 for the mean meridional flow near the surface. 

4.2 Solar cycle variations of meridional flow 

The meridional flow varies with the solar cycle. The resid- 
ual from the mean flows at each latitude is a ±5 m s^ 1 con- 
verging flow towards the active re gion belt t hat increases 



in strength with increasing activi ty (iGizon & Remp el 2008; 
Gonzalez Hernandez et al.ll2008h . 



4.4 Sunspots 

There is considerable observational evidence that sunspots 
have a large effect on waves. It is well documented that 
sunspots absorb up t o 50% of th e energy and shift the phase 
of incoming waves (IBraun et a l. 1988, 1992). Various de- 
pendencies on frequency and radial order were also found. 
These results were originally discovered using 
Fourier-Hankel methods. 

In the past few years, more evidence has eme rged that 

the w ave signal is do minated by surfac e effects (e.g. IBraun & Birch§ 
20081: [Cally 2009; G izon etalj [20091) . FuU waveform for- 
ward modelling has also revealed that a shallow sunspot 
model is sufficient to quantitatively explain the observed 
change in amplitu de and phase of the waves in 
caused by a sunspot ( Cameron et al.ll2010l) . Despite the ac- 
cumulating evidence that sunspots are large perturbations to 
the waves, wave-speed inversions below sunspots have tra- 
ditionally been done under simplifying assumptions (weak 
and non-magnetic perturbations). Such results should be 
treated with caution due to the inherent problems of attempt- 
ing a linear inversion involving non-linear perturbations to 
the waves. 

Additionally, the inversions do not account for source 
suppression or details of the measurement procedure (e.g. 
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Braun & Birch 2008, sho w that th e filtering is crucial). Fig- 



ure 19 from .Moradi et al.l (l2009ah directly compares the re- 
sults of various sunspot inversions. All of the methods, ex- 
cept the inversion for time-distance analysis, appear to be 
consistent with an increasing wave-speed in the top 2 Mm 
and small perturbations below this. However, this does not 
mean either result is correct and the cause(s) of the discon- 
gruity needs to be fully understood before a reliable inver- 
sion can be made. 

For a m ore comprehensive review of sunspot seismol- 
ogy see e.g. lGizon et a U ('2009| ), and for su nspot models in 
local helioseismologv .Moradi et al.l (l2009al) . 



4.5 Acoustic power and magnetic activity 

It has been known for some time that acoustic power is 
altered by the presence of magnetic field. Suppression of 
acoustic power is evident in strong magnetic fields consis- 
tent with reduced convective sources. This suppression af- 
fects the computed cross-covariances and thus the travel 
times. 

It has also been observed that higher frequency 
(5-6.5 mHz) wave power is enh anced in regions surround - 



ing active regions and sunspots ( Hindman & Bro wnlll99 
called acoustic halos. Recently, Schunker & Braun] (120 1( 
found that the halos tend to be stronger in regions of in- 
termediate strength horizontal magnetic field. From such 
a basic relationship between the acoustic power and mag- 
netic field vector, they are able to reproduce the acoustic 
power maps from the vector magnetic field to a high de- 
gree. They also found that the high frequency ridges of the 
power spectrum are shifted to higher wavenumber Several 
mechanisms have been suggested for generating the halo 
although none reproduce all of the observed characteris- 
tics of the halo. The mechanisms include mode conversion 



(iKhomenko & Colladosll2009l) . mode scatteri ng jHanasoge 



the physical mechanisms; //) with realistic numerical simu- 
lations of magnetoconvection where waves are naturally ex- 
cited by convection; and ///) with numerical simulations of 
linear wave propagation through a prescribed static back- 
ground. 

The models of i) were instrurnental for discerning the 
relevant phvsics. ' Spruit & Bogdanl ( ll992h proposed that the 
absorption could be caused by mode con version, and sub- 
sequent stu dies expl ored the physics (e.g . lCally & BogdanI 



I993II1997I) . In fact JCrouch et alJ (l2005h demonstrated that 



mode conversion could explain both the observed absorp- 
tion and phase s hifts and match t he quantitative measure- 
ments quite well dCallv et al.l2003h . These models suggested^ 
that the mode conversion occurs in the top 1 Mm below 
the surface, where the Alfve n speed is c omparable to the 
sound speed. More recently, ICallyi (12009 ) has shown that 
in umbrae thermal effects are likely to dominate the waves, 
whereas in the penumbra magnetic effect dominate. While 
these important results elucidate the mechanisms, now we 
require more realistic simulations where we can also test the 
analysis techniques and quantitatively measure the effects. 
Approach //) is fully physical, but very computer in- 
tensive. These kinds of simulations allow a study of the 
behaviour of granulation, and thus changes in the proper- 
ties of the waves, in the presence of magnetic fields. The 
most recent s i mulat ions including a sunspot can be seen in 
Rempel et al.l (120091). although limited helioseismic analysis 



has been published (IBirch et al.ll2009l) . 

The third approach is less computer intensive and allows 
for parametric studies of magnetic structures. A disadvan- 
tage is that these simulations require a stable background 
model. Several groups have developed numerical codes de- 
signed to model oscillations in the vicinity of sunspots using 
linearised MHD equations: 



2009 1 ) . trapping under magnetic field can opies (iKuridze et al.l 
2008h . and a change in convective scales ( Jacoutot et al.l2008b .| 



Chou et al.l (120091) analysed the spatial distribution of 



acoustic emission, absorption and suppression characteris- 
tics of two active regions and found that they all correlate 
spatially with the magnetic field, including the plage re- 
gions. 

(l2009h 



- The lAC MHD code (IKhomenko et all l2008h has been 
used for a variety of studies. Some of the relevant results 
include that seismic waves below a sunspot travel faster 
I than in the quiet-Sun and that high-frequency fast waves 
refracted by the magnetic field just above the surface are 
possi bly the cause of acoustic halos ( Khomenko & ColladosJ 



2009h . 



Zharkov et al. 



observe enhanced power in 
sunspot umbra at large line-of-sight angles, i.e. perpendicu- 
lar to the field. They model the propagation of waves 
through localised magnetic field and conclude that it is 
likely they are observing slow modes. 



5 Recent theoretical developments 

To understand the seismic signal observed on the surface of 
the Sun it is essential to model the propagation of waves 
through a magnetic field. There are three ways that this has 
been done: i) with simplified models aimed at understanding 



- The SLiM code ( ^Cameron et alJuOOTI) together wit h a 
carefu lly constructed stable background JSchunker e t al.J 
l2010ah and su nspot model is able to reproduce observa- 
tional results (iCameron et alJl2010l) . The sunspot model 
itself was also somew hat constrained by the comparison 
dCameron et al.ll2008h . The code has also been used to 
study the separate effects of the magnetic field and ther- 
mal effects of a su nspot. 

- The SPARC code (IHanasoge et al.ll2006l) has been used 
to study waves surrounding magnetic flux tubes, in par- 
ticular the acoust ic halo ([Ha nasoge 2009). 

- The SAC code dShelvag et al.ll2008l l2009i) was used to 
analyse the propagation and dispersion of waves in a 
non-uniform magnetic field. They also found that the 
largest contribution to travel time perturbations was due 
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to the thermal structure in weak to intermediate mag- 
netic field strengths, but that mode conversion in strong 
fields is important. 

All of these codes give qualitatively similar results. 

A careful analysis of the effects of stabilising a solar 
model against convection ha s been done using the SLiM 



code (ISchunker et al.l 1201 Obh . The stabilisation procedure 



changes the properties of th e model , and in turn alters the 
eigenmodes. ISchunker et al.l (1201 Obi) attempt to correct for 
changes to the eigenmodes by modifying the background 
model. Since the models are stable, a model of the convec- 
tive sources needs to be implemented as well as a model 
of the damping of the waves. All of these parameters need 
to work in concert to produce seismically solar-like simula- 
tions to provide realistic interpretations of the helioseismic 
observations. 



6 Discussion and future studies 

Further confirmation of the surface effects of strong mag- 
netic fields have been found in observations, casting fur- 
ther doubt on over-simplified inversions for wave-speed and 
flows below magnetic active regions and sunspots. The ef- 
fects of filtering, source suppression and mode conversion 
of the waves in strong surface magnetic fields have yet to 
be fully quantified for each of the analysis techniques. In 
turn, a concerted effort is underway to quantitatively sim- 
ulate the seismic signature of the Sun in the presence of 
magnetic fields to explore these effects. The simulations are 
still being fully analysed and in the near future a quantita- 
tive measure of the effects may be used to infer the correct 
subsurface properties of a sunspot. 

Future efforts in local helioseismology will focus on 
deeper inversions, particularly of the base of the convection 
zone, and steps tow ards this have be en taken using time- 
distance techniques dZhao et al.ll2009b . 

A challenging goal is to detect the emergence of mag- 
netic flux before it is observed a t the surface. Some progress 
has be en made towards this (e.g. lHartlep et al.l2010tlKom m et 
20091 I2OO8.) . but strong definitive signatures have not been 
observed. Future studies of the effects of magnetic flux emer-| 
gence will require systematic helioseismic analysis of many 
emerging active regions and realistic simulations. 

The HMI instrument onboard the recently launched 
SDO very soon provide increased resolution observations of 
the seismic Sun and provide accompanying vector magne- 
tograms on a regular basis. Thanks to the improved spatial 
resolution we anticipate being able to do local helioseismol- 
ogy closer to the limb (and poles). 
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